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Abstract 


This monograph provides an extensive list of formulas for airfoil 
polynomials. These polynomials provide convenient expansion func- 
tions for the description of the downwash and pressure distributions 
of linear theory for airfoils in both steady and unsteady subsonic flow. 


Symbols 

F 

f 

Jn 

K 

Pn 

pl/2 p-1/2 
1 n ? r n 

p(— 1/2, 1/2) p( 1/2, —1/2) 

r n : r n 

V 

T„ 

^m? tn 
U n 

Un 

W 

z 

£ 


hypergeometric function 
arbitrary function 
Bessel function of first kind 

kernel function of pressure-downwash integral equation 
Legendre polynomial 

associated Legendre functions of first kind 

Jacobi polynomials 

arbitrary polynomial; airfoil pressure 

associated Legendre functions of second kind 
Chebyshev polynomial of first kind 
airfoil polynomial of first kind 
Chebyshev polynomial of second kind 
airfoil polynomial of second kind 
airfoil downwash 
real variables 
complex variable 
real variable 


Introduction 

The so-called airfoil polynomials are used as expansion functions to compute the pressure on 
an airfoil in steady or unsteady subsonic flow. In kernel function aerodynamics (e.g., ref. 1), the 
unknown pressure distribution on the airfoil is expressed as the solution of the integral equation 


w 



p(OK(x-Od( 


in which the nondimensional variables are the known downwash w, the unknown lifting pressure 
p. and the kernel function K. The integral extends over the chord, which is normalized such 
that the interval [—1,1] extends from the leading edge to the trailing edge of the airfoil. The 
pressure has a square-root singularity at the leading edge and a square-root zero at the trailing 
edge. The kernel function involves an integral of Bessel functions and has both a first-order 
pole and logarithmic singularity at the point x = £ on the airfoil. The downwash, pressure, and 
kernel function are functions of the flow Mach number and the airfoil oscillation frequency in 
addition to the indicated dependence on the chordwise variables x or £. 



The airfoil polynomials are particularly convenient expansion functions for the pressure and 
downwash. For example, the pressure may be represented by the expansion 

= \l 2 Un ^ 

» 71=0 

in terms of the airfoil polynomials of the second kind u n (£). These polynomials possess many 
nice properties such as the interdigitation of the zeros of the polynomials of the second kind u n 
for the pressure and of the first kind t n for the downwash. Applications of these polynomials 
to subsonic lifting theory may be found in references 1-3. Also, they may be used to describe 
the chordwise pressure distributions in three-dimensional wing theory. In the three-dimensional 
theory, the spanwise pressure is typically expanded in a series of Chebyshev polynomials of the 
second kind U n {y). 

These polynomials are actually renormalized Jacobi polynomials as described, for example, 
in chapter 22 of reference 4. The present monograph roughly follows the organization of that 
chapter; however, it contains many additional formulas. Many of the formulas can be derived 
from the trigonometric definitions of the polynomials or from their expressions in terms of the 
more familiar Chebyshev polynomials. 

Formulas are included that relate the airfoil polynomials to many of the standard mathe- 
matical functions. Of particular use are the expressions for the singular integrals occurring in 
linear theory aerodynamics. Tables I and II summarize the coefficients for t n (x) and u n (x), 
respectively. 

Mathematical Properties 

1. Definitions 


cos ( n -f i ) arccos x 

tn(x) = - ^ 7 - J ^ 

cos ( ^ arccos x J 

sin [ n + \ ) arccos x 

»„(*) = — 1 ^ v — ^ 

sin f ^ arccos x J 


2. Orthogonality Relations 


tn(x)t m (x ) dx = 


I f 1 IH 

7T 7-1 V 1 — 

Ij.J = 


(m = n ) 
(m ^ n) 

(m = n ) 
(m 7^ n) 


2 



3. Recursion Formulas 


^n+l0*0 2x^77, (x) — 0 

2,XUji (x) + 'Utt,— i(x) — 0 
t n + 2 (x) - (4x 2 - 2 )t n (x) + t n _ 2 (x) = 0 
u n + 2 (x) ~ (4x 2 - 2)u n (x ) + u n - 2 (x) = 0 

4. Relations to Other Functions 
4-1* Chebyshev polynomials . 



-4*2. Jacobi polynomials . 

!,w=2 2 »(' 2 ;V I i><- 1A1/2) w 

M*)=2 2 "( 2 ;)“‘pi IA - 1/2 >M 

-4.5. Hypergeometric functions . 

<„(*) = + 

w n (x) = (2n + 1)F ^-n, n + 1; |; ~ 

<»(*) = (-l)”(2n + 1)F (-n, n + 1; |; 

«n(*) = (-l) n r(-n,n + l;|;i±^) 
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4-4> Legendre functions. 




u n (x) 




faW ~H)£(iTf f^w 

%(*) = ( n + I) V* ( I"Z"i) f*" 1 / 2 (x) 


5. Generating Functions 


oo 


XI 1 \ ^—~2 

1 — 2 xz + z L 

1 + 2? 

2xz + z 2 


72—0 

OO 


X M*)z n = — 


n = 0 


6. Rodrigues’ Formulas 


t n (x) = 


(— l) n 2 n n! /I -a: 


( 2 n)! V 1 + x dx n 


i^(l-x 2 )” 

1 — X 


Unix) = 


(— l) n 2 n n! fl + x d n 


( 2 77,)! V 1 — x dx n 


1 — X 

l + x 


2\n 


(i - **) 


7. Interrelations 


U n {x ) = (-1 ) n t n i~x) 

(1 + x)t n ix ) = XUnix) — U n -lix ) 

(1 + x)t n ix) = -XUnix) + U n+ lix) 

(1 - x)u n ix) = -Xtnix) + tn-lix) 

(1 - x)u n ix) = Xtnix) - t n+ lix) 
Unix)t n -lix) = t n ix)u n -lix) + 2 
tnix) tn — 1 (^) — Unix) Uji — i(ar) 

(1 + X) [-2nOO + *n-l(z)] = (1 - X) [Unix) + U n -\ix)\ 
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8. Explicit Expansions 


*»<*) = 2_ “ E(-i)* ( 2 ” 2 £ 1 ) a + - *)‘ 

k = 0 V J 

“»(*) = 2_ ” E(-l) i f « + I 1 ) (1 + *)"-*(! - *)* 
&=0 V 7 

^ = 2-"E( l( „_Vj)^ ) 

■• = r * ( L( n-% /2J )“‘W 

9. Expressions for Negative Degree 

^ — 72— 1(*^) = ^72 (*^) 

= - u n (x ) 

10. Differential Equations 



11. Hilbert and Related Transforms 
11.1. Hilbert transforms . 


I f 1 f l + x t n( x ) 

n J-i V 1 — £ x — y 

1 f 1 l l ~ xu n(x) 

n 7-i V 1 + x x - y 


dx = u„(y) 
da; = -t n (y) 
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1 f 1 /I — a; iS(a;) , 1 — y . . . . 

-<*« = - — - t n (y)u n (y ) - 


7 T J_1 V 1 + £ £ — 2/ 


l + y 



1 -f- 2/ 

dx - -- t n (y)u n (y ) + 

1 - 2 / 


1 — 71 
eta = / n (y) In - : - + 4 


(n ± 1) 

(n 7^ 1) 

2(2n + l) 

1 + 2 / 

2(2n + l) 


/ 


1 t/ ra (x) 
'-1 * - 2/ 


1 + 2 / 

da; = ■u„(t/)lnf- r ^ + 4 


1 - 2 / 

^n— fc( 2/) 

it 


£= 1 , 3 ,... 

1 - y , . NT' u n-k(y ) 



1 + a; / n (a;)p(a;) 


7T 7_! V 1 - a; x-y 

1 / l-a;u n (a;)p(a;) _ 


1 



1 + 2 / 

da; = u n (y)p(2/) 
dx = - t n (y)p(y ) 


k= 1 , 3 ,... 


k 


7T 7_! V 1 + * x - y 

where p(x) is any polynomial of degree n or less. For complex z, 

1 f 1 /I + x t n (x) 


ft J - 1 V 1 — X X — z 


1 


7 r 7-i V i + x x — z 



1 - x ti„(a:) 


dx = u n (z) 


dx = -/ n (z) + 


2+ 1 , s 

- r^nl^) 
z — l 


z — 1 

* + l 


u ra (z) 


j. f 1 1 

* J - 1 Vl - a ; 2 


dx = 7=^ . 

^(x — z) \/^ — 1 


11 , 2 . Finite part transforms. 


1 f 1 / 1 + X tn{x) J -(2 n + l)/n( 2 /) + «n( 2 /) / / , 

? L V Tziorr^i * = 2 (T^) = “»« 

I f 1 / l-^ ^0*0 ^ _ ~( 2n + 1 ) tf n(2/) + ^n(y) _ _ , , 

7T y_i V 1 + « (x - 2/) 2 * ~ 2(1 + y) ~~ 
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I [ 1 / l~g(l + g)*n(g) , _ — (2n + l)< w (y) - u n (y) 

-k J_i V 1 + x (x - y) 2 X ~ 2 

1 f 1 ll + x(l-x)u n (x) _ ~(2n + l)u n (y) - t n (y) 
W-iV 1-x (x-y) 2 2 


11.3. Logarithmic transforms . 



1 + £ 


7T J_i V 1 “ S 


i n (x)ln|o? — 2/| dz = 


«n-l(y)-«n(y) «n(y) - «n+l(y) 

2n 2(n + 1) 


- In 2 - y 
_ r n+ i(y) T„(y) _ f M 

iT+I — - jMv)iv 



1 — X 


1 T 7_1 V 1 + X 


u n (x ) In | a; — y| da? = < 


f t n -i(y) + <n(sO <»(?) + Wi(y) 


2n 

l - In 2 + y 


2(n + 1) 


= 7V^)_TUy) = r tn{)d 

n + 1 n J 


12. Derivatives 



(1 ■ c ) [^n— 1(*) H“ ^n( a: )] — n [^n— lC 2 -) ^n(*)] 


(n ^ 0) 

0 = 0) 

(n^O) 
(n = 0) 
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n—1 


4iP ) = - E^ 1 )” k [ n ( n + *) _ k ( k + ^WP) 

k=0 


n—1 


* 40 ) = E ^ 71 + !) - k ( k + !)]^ 0 ) 

>b=0 


n—1 

40 ) + 4 -iO) = 2n E OP) 

k = 0 


n—1 


«nO)-«n-lO) = - 2n E(“ 1 )" *“*(*) 

fc =0 


13. Indefinite Integrals 


/ ^737^0) = 2n(n+"lj '^ B ( g ) ~ ^ 2n + 

—01 — 


= < 


Un(x) U n - lO) 


l. — arccos a; — \/l — a; 


n + 1 n 


p > 0) 

(n = 0) 


Y X yj\ — X ‘2 

u n (x) dx = , ^ [-MnO) + ( 2n + l)^n(*)] 


1 + X 


2 n(n + 1) 


= < 


o„0) U n —i 0) " 


y/l-x 2 , 

n + 1 n 

l — arccos a; + Vl — a: 2 


p > 0) 

(n = 0) 



*n-lQ)+*np) <np) + *n+lP) 
2n + 20 + 1) 

U n - lQ) ~ u n {x) Un Q) - U n +lP) 

2n + 20+1) 


4»0) , ^~n+lQ) 
n n + 1 

7nQ) _ ^n+lp) 
n n + 1 


| X — - — da; = arccos a; — \ 1 In 

1 — a: x — y \ 1 — y 



-In 


1 - * /I - y 


1 + x 


1 + 2 / 
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14. Cross Expansions 


n 

u n(x) = 2 ^ ^ t m (x) — t n (x) 
m= 0 

n 

tn(x) = 2 ]T (-1 ) n - m Um(x) - u n (x) 

m = 0 


15. Multiplication Formulas 

In all cases n > m. 

2t n {x)Um(z) ^n+m(^) “1“ ^n+m(^) V? n—m{x ) + t n—m (*) 

2 U n {x)t m {x) = U n +m(x ) + t n + m {x ) + U n - m ix ) — t n -mix ) 
2t n (x)t m ( < x) = t n _f_ m (x) + — m ( x ) + ^ ^ [ w n+m(^) + ^n-m(^)] 

2(1 + ^)tn(^)tm(^) = £77+777(2:) + £77+77741(2') 4 £77-777(2:) + £ n _ m _i 

2w n (x) / a m (x) = ix n -f m (x) 4 U n - m (x) 4- —[—t n +m(x) 4 t n - m (x)\ 

2(1 — x)u n (x)u m (x) = ^n+m(^) ~ ^n+m +1 (^) + U n - m (x) — 7X n _ m _i(x) 


16. Christ offel-Darboux Formulas 


n 

^ ^ u m(x)u m {y ) 

777-— 0 


u n +i(x)u n (y) - Un(x)u n +i(y) 
2(x - y) 


n 

y! tm(x)t m (y) 

771=0 


^n+1 (jQ^n(y) ^77(^)ty l -{-i(y) 

2(x - y) 


n 

^ ^ 777 (^)^ 777 ( 2 /) 

777=0 


2 4 * n +l(^)^n(y) ~ ^n(x)^ n+ i(y) 
2(x - y) 


y! ^777 ( X ) — 9 [' U 77+l( :r )' U n(^) — tt 7l(^)Wr7+l(^)] 
777=0 


^ 1 

^ V ^777 (*^0 “ 2 [^W+l( a ')^w(^) ~ ^77(^)^7741 (^)] 
777=0 
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tm(x)u m (x) = ^[ 4 + l(^)u n (x)-t f n (x)u n+1 (x)] 


x-y 


2 [ U n+l( x )tn(x) U n (x)t n -\- l(#)] 

oo 

= 'y ^ tm(x)u m (y) 


17. Special Values 


tn(~~ 1) — (— l) n (2n + 1) 

t n \y) = ±1 = cos — sin — 

2 2 

*n(l) = 1 

M-i) = (-i) n 

/«\ , T17T .7171 

U n (0) = ±1 = cos 1- sin — 

2 2 

u n (l) = 2n + 1 


18. Zeros 


t n (xi ) = 0 where X{ = 


u n (xi) = 0 where X{ = —> 


-7T (* = 1,2,..., n) 


2i — 1 


jtt (* = 1,2, ... ,n) 


19. Miscellaneous Definite Integrals 


1 ( 1 (m = n = 0orm = n + l) 

^ J y/l - x 2 t n (x)it m (x) dx = < -1 (m = n - 1) 

1 0 (otherwise) 

1 j 1 -j- x 

- J Y x u n(x) u m(x) dx = 4min(m, n) + 2 - <$ mn 

d:c = (- 1 ) n_m [ 4miii (”*,«) + 2 - <WI 
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n — m 


/ f n (a:)« n (x) dx = 

J - 1 2n + 1 

J u n (x) dx = 


n + m + l 

2 


n + 1 


l + (~l) n + l-(~l) n 


n + 1 



1-x. 


v i + x 1(I - y )n dx - D-d-*-* (*)( fc / 2 ) s'” -1 


k=0 


1 Z" 1 1 — xy + Vl — x 2 \/l — y 2 

V 1 -y 2 


1 y 1 ^ 1 — xy + VI — x 
77 J - 1 1 — xy — Vl — x 


t n (x) dx = 


r 1-y 


Uniy) ±Un+]Xy ) «n-l(y) + «n(y) 


1 + y L ra + 1 

2(1 +y)y] 


n 


/ 1-y 

l + y 


20. Relations Involving T n (x) 

tn(x) ”f~ tfi — i(x) — Unix) Uji — j(x) — 2Txi{x') 

+ /-"v _ Tn{x) + T n +\(x) 

‘nw — z ; 


(n # 0) 


1 + x 


Unix) = — Tn + l(x) 
1 — X 


^n(^)^n(^) 1 (^)^n— 1 (^) — 2^2 n(*^) 


(«^o) 



<n(x) = r„(x) - (1 - x)£/ n _i(x) 
«n(x) = r n (x) + (1 + x)U n - i(x) 


(n > 0) 

( n > 0 ) 

(n = 0) 
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21. Relations Involving U n (x) 


tn(x) = U n (x) - U n -i(x) 
u n(x) Unix) + U n -lix) 
tnix) + Unix) t n -\-\ix) + U n -\-\ix) — 2U n ix) 

tnix)u n ix) = U 2n ix) 



t n ix) = Tnix) - (1 - x)U n - l(x) 
Unix) = T^x) + (1 + x)U n -lix) 

tnix)~t n+ i(x) U n ix) + Un+lix) nTT , s 

= TTx = 2tr »W 

22. Integrals Involving Legendre Polynomials 


tnix) - 


Unix) = 


2n + l f x Pniy) 


2\/l + x J - 1 y/x-y 

2n + l f 1 P w (y) 
2\/l Jx v ll ~x 

tniy) 


Pnix) = - f 1 7f =L 

^ da; V 1 ~ 2/Vy ~ ® 

p»(x) = - r un(y) 

n j-i 


-i V 1 + vV x - y 


dy 

dy 

dy 

dy 


23. Fourier Series 


f( x ) = Yl tnix) 1 / ^^7^n(y)/(y) dy 

n=0 ' /_1 V y 

1 /*1 /J 

/(*) = 2 Un ( x ) ~ / V T+T u n(y)fiy) dy 

n= 0 • /_1 V y 

h _l x 1 /*l 

/(*) = V 1 — / tniy)fiy)dy 

' 1 ~ X nZo n J~ 1 

/ "][ — x 1 rl 

/( x ) = v 777 2_/ “»(*)“ / «n(y)/(y) dy 

Vi + x n =o 71 - 7 - 1 
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24. Finite Fourier Transforms 


l j\ e " y Vr^ n(x) dx = injn{y) + * n+lj -+^) 

£ /i e “ y vXXf Un(rc) dx = * n Jn(y) " * n+1 J -+ i(y) 

- / 1 e* x?/ \/l - x 2 t n (x) dx = + Z>+l) Jn+l(y) 

7 T 7 -i y y 

— [ e lxy Vl - X 2 u n {x) dx = — J n (y ) + + ^ J n+1 (y) 

v J-i y y 

Numerical Methods 
1. Gaussian Quadrature 


7r Xi Xl - dX ~ 1 + 


2n + 1 f , 
2=1 


where x 7 = — 


2i 

COS 7T 

2n -b 1 



, ./ T — — f(x)dx= - - - y^(l — Xi)f(xi) 

7T 7_! V 1 + x 2n + l^ v l> 

2=1 


, 2i - 1 

where x; = — cos 7 r 

1 2n + l 


2. Hermite Quadrature 


\ !\'Ir£ ,(x)ix 


2n -j- 1 


[(1 + Xi)f(xi ) - 2(1 - x ;)/'(x;)] where x* = - cos - n 

i — 1 


\ S\'lWx i{x)ix 


^ (1 - Xi)f(xi) + 2(1 - x?)/'(x;) where x t = -cos 


2n + l ^ 


i f i jm 

TT 7-1 V 1 + x 


U m (x)f(x) dx 


2n + 1 , 

2 = 1 


y> + Xi) [(2m + l)i TO (x;)/(x;) - 2(1 - Xi)um(xi)f'{xi)] where x z = 


2 i 

— cos -7 r 

2n + 1 
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, 2 * + 1 

where x; = — cos n 

2n+ 1 


- f J\^-t m {x)f{x)dx 
7T 7-1 V 1 — X 

2 n 

-stiB 1 - Xi) [(2m + l)ttm(xi)/(x,) + 2(1 + x i )i TO (x i )/ , (x i )] 

Z = 1 


Concluding Remarks 

An extensive list of more than 150 formulas for the so-called airfoil polynomials has been given. 
These orthogonal polynomials are convenient for describing the chordwise pressure distribution 
in linear theory aerodynamics. Of particular use are the expressions for the singular integrals 
occurring in the linear theory. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
September 30, 1994 
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Table I. Coefficients for t n (x) 

" t n (x ) — c m% m and x n — b^ 1 Y^m=o dmtm(x)', values in table correspond to" 

coefficients as follows: Cm are given horizontally below the diagonally stepped line, dm 
are given vertically above the diagonally stepped line, and b n as indicated for a given 
x n ; e.g., tj(s) = 1 + 4s — 12s 2 — 8s 3 + 16s 4 and s 4 = yjj(6to + 4ii + 4f2 + + 14 ) 



n 

X 

X 1 

s 2 

2 

X 

X 4 

s 5 

X 6 

s 7 

X 

1 0 

i 

1 

2 

3 

6 

10 

20 

35 

70 

i 




*1 

-1 

1 

1 

3 

4 

10 

15 

35 

56 


2 




1 2 

-1 

-2 

i 

1 

4 

5 

15 

21 

56 



4 




■ 

1 

-4 

-4 

i 

1 

5 

6 

21 

28 




8 




D 

1 

4 

-12 

-8 

i 

1 

6 

7 

28 



16 




h 

-1 

6 

12 

-32 

-16 

i 

1 

7 

8 



32 



^6 

-1 

-6 

24 

32 

-80 

-32 

i 

1 

8 



64 


*7 

1 

-8 

-24 


80 

-192 

-64 

i 

1 




h 

1 

8 

-40 

-80 



-448 


1 


256 

bn 

i 

2 

4 

8 

16 

32 

64 

128 

256 


15 























































Table II. Coefficients for u n (x) 

~ u n {%) = Y^m=Q c mX m and x n = b^ 1 X)m=0 d m Um(xy, values in table correspond to' 
coefficients as follows: c m are given horizontally below the diagonally stepped line, d rn 
are given vertically above the diagonally stepped line, and b n as indicated for a given 
x n ; e.g., u^x) = 1 — 4x — 12a; 2 + 8a; 3 + 16a; 4 and a; 4 = -^(6uo — 4ui + — u% + U4) . 



X® 

x 1 

X ^ 

a; 3 

X 4 

a; 5 

a; 6 

X 7 

a; 3 


1 

-1 

2 

-3 

6 

-10 

20 

-35 

70 




Ul 

1 

1 

-1 

3 

-4 

10 

-15 

35 

-56 

2 



U2 

-1 

2 

1 

-1 

4 

-5 

15 


56 


4 



u z 

-1 

-4 

4 

1 

-1 

5 

-6 

21 

-28 



8 



U4 

1 

-4 

-12 

8 

1 

-1 

6 

-7 

28 


16 





1 

6 

-12 

-32 
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